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20.  ABSTRACT  (Continued) 

'4^The  application  of  the  CAPRI  code  to  numerically  analyze  the  response  of  a 
magnesium  nose  capped  cylindrical  steel  projectile  striking  sandstone  is 
demonstrated  and  discussed.  Detailed  numerical  results  depicting  the  config- 
uration of  the  projectile-sandstone  system,  and  contours  of  stress  in  the 
projectile,  at  various  times  after  impact  are  given.  Also,  time  histories  of 
acceleration,  velocity,  displacement  and  stress  at  selected  material  locations 
are  shown. 

The  CAPRI  code  represents  an  initial  attempt  to  develop  an  Eulerian  numerical 
code  suitable  for  numerical  analysis  of  such  a class  of  impact  problems.  It 
is  suggested,  however,  that  further  refinements  be  made  in  several  features 
of  the  numerical  scheme  before  the  code  is  made  generally  available  for 
projectile  design  studies. 
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PREFACE 


It  has  been  suggested  that  one  possible  way  of  cushion- 
ing the  impact  of  an  earth  penetrating  weapon  (EPW)  striking 
rock  at  high  velocity  may  be  to  add  a sacrificial  nose  cap  to 
the  EPW.  Such  a nose  cap  would  be  designed  to  undergo  large 
plastic  flow  and  fracture  during  the  impact  process  and,  in  so 
doing,  dissipate  a significant  portion  of  the  initial  kinetic 
energy. 

The  investigation  reported  herein  was  part  of  a larger 
effort  sponsored  by  the  Defense  Nuclear  Agency,  Strategic 
Structure  Division,  to  examine  the  feasibility  of  the  sacrific- 
ial nose  cap  concept.  The  objective  of  Systems,  Science  and 
3 

Software's  (S  ) participation  in  the  program  was  to  develop 
a two-dimensional  Eulerian  numerical  method  suitable  for  analy- 
zing EPW-rock  impacts  which  could  be  used  for  parametric 

studies  to  perfect  the  design  of  sacrificial  nose  caps.  The 

3 

work  performed  by  S toward  this  goal  is  described  herein. 

The  Principal  Investigator  for  S^  was  Dr.  H.  E.  Read,  and  the 
Project  Officer  for  the  DNA  was  Lt.  Col.  D.  Spangler. 
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throughout.  Thanks  are  also  extended  to  Mr.  M.  H.  Wagner, 
California  Research  and  Technology  (CRT) , for  helpful  conversa- 
tions pertaining  to  the  CRT  sandstone  model  and  for  performing 
a Lagrangian  stress  propagation  calculation  using  this  model  to 
aid  us  in  checking  our  algorithm  for  incorporating  the  model 
into  our  Eulerian  code. 
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1 . INTRODUCTION 


1 . 1 BACKGROUND 

Considerable  interest  has  developed  within  the  defense 
community  during  recent  years  in  weapons  that  are  capable  of 
penetrating  some  distance  into  the  earth  before  detonating. 

Such  a weapon  function  is  particularly  attractive  because  it 
(a)  reduces  the  collateral  damage  produced,  and  (b)  enhances 
the  direct-induced  ground  shock  developed.  To  be  generally 
effective,  earth  penetrating  weapons  (EPW)  must  be  capable  of 
penetrating  various  geologic  targets,  ranging  from  soil  to  rock 
without  producing  critical  damage  to  the  payload.  To  design 
weapons  with  such  a mission  capability  requires  knowledge  of 
the  intensity  of  the  dynamic  loads  experienced  by  the  reentry 
vehicle  and  its  internal  components  during  the  penetration 
process . 

The  major  effort  to  date  to  understand  projectile 
penetration  reponse  has  been  directed  toward  soil  targets.^  ^ 
In  this  case,  the  hard  (usually  steel)  penetrators  undergo 
relatively  small  deformation  during  penetration  into  soil.  An 
extensive  data  base  on  projectile  penetration  has  been  accumu- 
lated over  the  years  from  numerous  field  tests  involving  both 
large  and  small  scale  penetrators  and  a variety  of  soil  targets 

The  penetration  of  projectiles  into  hard  rock  is  also 
of  defense  interest  but  has  received  less  attention  — both 
experimentally  and  numerically  — than  soil  penetration. 

In  an  effort  to  reduce  the  intensity  of  the  stress  waves  that 
are  propagated  back  into  the  main  body  of  such  a penetrator 
after  impact  on  rock,  several  modifications  to  the  basic  pene- 
trator design  have  been  proposed.  One  of  these  involves 
attaching  a nose  cap  of  different  material  to  the  penetrator. 
The  notion  behind  this  concept  is  that  the  nose  cap  material 
can  be  designed  to  undergo  large  plastic  flow  and/or  fracture 
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during  the  penetration  process,  and  in  so  doing,  dissipate 
sufficient  kinetic  energy  tc  cushion  the  impact  process. 

In  the  case  of  projectiles  which  are  sufficiently  hard 

that  they  experience  very  little  deformation  on  impact  with 

( 8 ) 

rock,  the  Lagrangian  numerical  method'  originally  developed 
for  projectile-soil  penetration  appears  to  be  adequate.  How- 
ever, if  on  impact  with  rock  a projectile  experiences  large 
deformation,  the  Lagrangian  numerical  approach  described  in 
Reference  8 is  inadequate.  For  such  cases,  a numerical  method 
is  needed  that  can  cope  with  large  plastic  deformations  accur- 
ately, and  can  also  treat  the  effect  of  sliding  friction, 
which  develops  between  the  penetrator  and  the  rock,  on  penetra- 
tor  deceleration. 

Unfortunately,  none  of  the  existing  Lagrangian  or  Eulerian 
numerical  methods  are  suitable  for  this  class  of  impact  prob- 
lems. While  Lagrangian  methods  are  capable  of  treating  mul- 
tiple materials  and  frictional  sliding  quite  accurately,  they 
are  unsuitable  for  problems  involving  large  plastic  flow;  al- 
though Lagrangian  methods  are  often  applied  by  the  defense 
community  to  problems  involving  large  plastic  deformations, 
the  calculational  process  usually  requires  extensive  rezoning, 
and  little  is  known  about  the  credibility  (accuracy)  of  the 
results  so  obtained. 

The  basic  strength  of  the  Eulerian  method  lies  in  its 
ability  to  treat  large  deformations  straightforwardly.  On 
the  other  hand,  the  Eulerian  method  is  not  particularly  well- 
suited  for  handling  explicit  material  interfaces,  especially 
those  defining  free  surfaces  and  sliding  interfaces.  The 
Eulerian  approach  is  also  susceptible  to  spurious  numerical 
diffusion,  which  arises  from  motion  of  material  through  a fixed 
Eulerian  grid.  In  Lagrangian  schemes,  on  the  other  hand, 
the  grid  is  attached  to  and  moves  with,  the  material  and  con- 
sequently there  is  no  mass  transport  across  cell  boundaries. 
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The  mass  transport  procedure  typical  of  most  Eulerian  schemes 
requires  a mixing  between  the  advected  mass  and  that  which 
remains  in  a zone;  this  mixing  is  responsible  for  the  numerical 
diffusion,  but  dees  tend  to  stabilize  the  numerical  scheme. 

While  the  Lagrangian  and  Eulerian  methods  each  have 
their  own  characteristic  limitations  and  virtues,  we  decided 
to  adopt  the  Eulerian  approach  in  the  present  effort  as  the 
basis  on  which  we  would  attempt  to  develop  a numerical  method 
suitable  for  the  analyses  of  penetrator-rock  impact.  Because 
of  the  severe  distortion  expected  to  occur  in  the  class  of 
problems  of  interest  to  the  present  study,  it  was  felt  that 
the  Eulerian  approach  would  be  preferable  if  the  errors  due 
to  numerical  diffusion  did  not  significantly  degrade  the 
numerical  solution.  Furthermore,  an  improved  mass  transport 
algorithm  would  be  needed  to  reduce  numerical  noise  generated 
as  material  interfaces  cross  cell  boundaries.  It  was  recog- 
nized at  the  outset  that,  to  accomplish  this  objective  would 
require  a number  of  significant  advancements  in  Eulerian  code 
technology,  and  the  present  study  was  undertaken  with  this 
risk  in  mind, 

1.2  SCOPE  OF  EFFORT 

The  scope  of  the  effort  described  herein  was  limited  to 
(a)  the  development  of  a two-dimensional  numerical  method  suit- 
able for  analyzing  the  impact  of  projectiles  with  rock  targets, 
and  (b)  solution  of  a problem  involving  the  normal  impact  of 
a solid  projectile  having  a soft  nose  cap  with  hard  rock  to 
demonstrate  the  feasibility  and  applicability  of  the  new  num- 
erical method.  Attention  was  confined  to  solid  projectiles 
since  it  was  recognized  that  further  advancements  in  numerical 
methods,  beyond  the  scope  of  the  present  effort,  would  be  re- 
quired to  produce  a procedure  that  would  be  satisfactory  for 
analyzing  the  impact  reponse  of  soft  nose  capped  penetrators 
having  thin-walled  main  bodies. 
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1.3 


SUMMARY  OF  RESULTS 


The  main  objectives  of  this  study  were  accomplished, 
namely  (a)  a new  Eulerian  code  CAPRI*  was  developed  which  is 
specifically  tailored  to  treat  problems  involving  the  normal 
impact  of  deformable  projectiles  with  rock,  and  (b)  the  appli- 
cation of  CAPRI  to  a projectile-rock  impact  problem  was  demon- 
strated . 

The  CAPRI  code  is  a two-dimensional  Eulerian  code  devel- 
oped from  the  existing  HELP  code  ^ framework.  It  features 
I the  following  new  capabilities  which  are  evidently  important 

for  the  class  of  impact  problems  of  interest  to  the  present 
investigation : 

• Accounts  for  material  strength  in  mixed- 
material  computational  cells. 

• Treats  frictional  sliding  between  material 
surfaces  in  contact. 

• Satisfies  the  correct  boundary  conditions  on 
stress  and  velocity  at  interfaces  between 
materials  in  contact  by  equilibrating  the 
normal  components  of  stress  and  velocity, 

• Contains  an  improved  algorithm  for  trans- 
porting mass  and  other  related  material 
properties  through  the  Eulerian  grid. 

• Accommodates  complex  geologic  material  models, 
including  those  having  a number  of  history- 
dependent  parameters  (see  Appendix  A for 
further  discussion  on  this  point) . 

The  major  portion  of  study  was  devoted  to  the  formulation, 
development  and  debugging  of  the  intricate  numerical  algorithms 
required  to  incorporate  the  new  features  listed  above. 

To  examine  and  demonstrate  the  ability  of  the  new  code 
CAPRI  to  perform  numerical  analysis  of  deformable  projectile- 
rock  impacts,  it  was  applied  to  a problem  involving  a solid 


*CAPRI  is  an  acronym  for  Code  for  Analysis  of  Projectile-Rock  j 

Impacts. 


cylindrical  projectile,  having  a soft  magnesium  nose  cap, 
impacting  sandstone  at  normal  incidence.  A numerical  solution, 
which  appears  to  be  acceptable,  was  obtained  out  to  a time 
after  impact  (27  usee)  at  which  a substantial  part  of  the  pene- 
tration process  had  occurred.  In  the  process  of  performing 
this  calculation,  however,  a number  of  difficulties  with  the 
numerical  scheme  arose  which  made  it  necessary  for  us  to  discon- 
tinue the  calculations.  Each  time  this  occurred,  the  source  of 
the  difficulty  was  uncovered.  An  appropriate  modification  was 
made  in  the  flow  field  at  that  time*,  and  the  calculation  resumed. 

The  calculation  was  not  continued  beyond  27  usee  because, 
by  this  time,  the  soft  magnesium  nose  cap  material  had  undergone 
such  severe  plastic  flow  that  its  thickness,  in  some  regions, 
became  less  than  the  Eulerian  zone  size;  this  resulted  in  three 
materials  being  present  in  one  zone  and,  in  the  present  state 
of  development,  CAPRI  (as  well  as  most  Eulerian  codes)  is 
unable  to  adequately  cope  with  such  a situation.  However,  the 
numerical  calculation  could  have  been  rather  straightforwardly 
continued  by  rezoning  the  Eulerian  grid  with  smaller  zones; 
this  was  not  done. 

The  work  reported  here  clearly  represents  only  a first 
step  toward  developing  a numerical  method  specifically  designed 
for  analysis  of  projectile-rock  impacts.  To  adequately  incor- 
porate some  of  the  physical  effects  felt  to  be  important  in 
such  impact  problems  into  the  Eulerian  framework  possibly  re- 
quires extending  the  Eulerian  method  to  the  limits  of  its  cap- 
ability. Clearly,  there  is  a need  for  further  investigation  of, 
and  improvements  in,  a number  of  features  of  the  CAPRI  code 
before  it  is  made  available  for  general  use  in  penetrator  design. 

*Each  time  a modification  in  the  flow  field  was  made,  it  was 
confined  to  a single  cell.  Furthermore,  the  cells  in  which 
such  modif ications  were  made  were  located  in  regions  where 
their  effects  probably  had  negligible  influence  on  the  over- 
all deformation  and  penetration  process. 


2.  NEW  FEATURES  OF  THE  CAPRI  CODE 


The  major  new  capabilities  which  are  available  in  the 
CAPRI  code  are  described  in  this  section;  these  include:  (a) 
an  approach  for  treating  frictional  sliding,  (b)  a procedure 
for  incorporating  material  strength  in  mixed-material  zones, 
and  (c)  an  improved  algorithm  for  transport  of  mass  and  other 
related  material  properties. 

2.1  FRICTIONAL  SLIDING 

Eulerian  numerical  procedures  have  been  developed  and 

incorporated  into  codes  for  treating  frictionless  sliding  be- 

( 9 ) 

tween  materials  m contact.  In  many  problems  of  practical 

interest  which  involve  the  sliding  between  two  metals  in  contact, 
the  heat  generated  by  the  sliding  process  may  produce  melting 
at,  and  in  the  vicinity  of,  the  material  interface.  In  many 
such  cases,  it  can  be  reasonably  assumed  that  the  sliding  pro- 
cess is  frictionless  and  this  simplification  leads,  in  many 

instances,  to  predictions  that  exhibit  reasonable  correlation 

(9) 

with  observations. 

For  the  class  of  problems  of  interest  to  the  present 
study,  however,  it  was  felt  that  the  frictional  forces  due  to 
sliding  may  have  a significant  effect  on  nose  cap  deformation 
and  overall  penetration  performance.  On  this  basis,  a numerical 
procedure  for  treating  frictional  sliding  was  developed,  in- 
corporated into  the  CAPRI  code,  and  utilized  in  the  numerical 
study  discussed  later  in  Section  3. 

In  the  sequel,  the  term  "slip  cell"  is  used  to  denote  an 
Eulerian  computational  cell  which,  at  the  time  of  interest,  con- 
tains two  materials  that  may  slip  relative  to  each  other  along 
their  common  interface.  The  sliding  materials  may  be  either 
similar  or  dissimilar.  In  the  sections  which  follow,  the  pro- 
cedure developed  for  treating  sliding  friction  in  the  CAPRI 


code  is  described.  Parenthetically,  we  note  that  the  Eulerian 
numerical  scheme  on  which  CAPRI  is  based  consists  of  three  basic 
phases,  namely,  (1)  a Lagrangian  phase  in  which  the  governing 
equations  are  treated  by  disregarding  the  effect  of  mass  flow, 

(2)  a flow  phase  in  which  the  material  is  transported  through 
the  fixed  Eulerian  grid,  and  (3)  a final  phase  in  which  the 
material  components  in  a cell  are  homogenized.  The  major 
revisions  to  the  current  Eulerian  methodology  were  made  in  the 
Lagrangian  and  Mass  Transport  phases,  and  the  details  of  how 
this  was  accomplished  are  given  below. 

2.1.1  Equat ions  for  Lagrangian  Phase 

In  the  CAPRI  numerical  scheme,  the  material  interface 
is  defined  by  a string  of  (passive)  tracer  particles  which  are 
attached  to,  and  move  with,  the  material  as  it  travels  through 
the  fixed  grid  during  the  deformation  process  (see  Figure  1). 

To  facilitate  the  treatment  of  sliding,  the  interface  is  con- 
sidered to  be  defined  by  a straight  line  connecting  the  points 
where  the  tracer  string  intersects  the  boundaries  of  the  slip 
cell,  i.e.,  line  ab  in  Figure  2.  The  normal  and  tangential 
directions  are  defined  in  terms  of  the  angle  0 which  the  line 
ab  makes  with  the  radial  axis. 

Figure  3 depicts  the  stresses  acting  on  each  of  the  two 

material  elements  occupying  a typical  slip  cell.  Across  the 

common  material  interface  the  tangential  component  of  velocity 

, is  in  general,  discontinuous,  while  the  following  continuity 

conditions  hold  for  the  normal  stress  o , the  tangential  shear 

n 

stress,  t,  and  the  normal  component  of  velocity  ; 


Figure  1.  Mixed  Eulerian  ceil 
with  slip  interface. 


Figure  2.  Idealization  of  slip 
interface  employed  in 
sliding  algorithm. 


Figure  3.  Stress  components  acting  upon 

material  elements  in  a slip  cell. 
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Here,  and  in  that  which  follows,  the  subscript  1 refers  to 
material  1 and  subscript  2 to  material  2.  Upon  applying  the 
law  of  motion  over  a time  interval  At  to  each  of  the  material 
elements,  the  following  system  of  linear  equations  result  which 
govern  the  motion  of  the  material  elements  during  the  Lagrangian 
phase  of  the  numerical  scheme: 


/M 

sin 

OA. 

\a 

- I 

/.At 

cos  0A. 

I.  - /u? 

+ 

— f"  \ 

(2) 

K 

l 

I n 

! 

iml 

l 

1 

' l 1 

ml  1) 

f- 

cos 

OA. 

\o 

- / 

f At 

sin  OA- 

V ’ (v? 

+ 

(3) 

\mi 

1 

I n 

l 

I \ 1 

ml  2 / 

/At 

sin 

OA. 

\ 0 

+ / 

At 

cos  OA. 

•4 

II 

c 

) 3 

+ 

At  Fn  V 

(4) 

V2 

l 

/ 

1 n 

\ 

m2 

l 

l \ z 

m2  3 / 

/At 

cos 

OA. 

\o 

+ / 

' At 

sin  OA. 

\t  - /v" 

+ 

At  Fn  v 

(5) 

\m2 

1 

1 n 

\ 

m2 

l 

/ \ 2 

m2  4 ) 

ml 

and 

m 2 1 

denote 

the 

masses 

of  the  material  elements 

occupying  the  slip  cell  at  time  t , and  A^  is  the  surface 
area  of  the  interface  between  the  materials  at  tn.  The  tilde 
is  placed  above  the  velocity  symbols  to  denote  their  values  at 
the  end  of  the  I.aqrangian  phase.  Finally,  the  forces  F1^ 
are  defined  by  the  following  equations: 


^rr’lR  A1R  " (orr)lL  A1L  + (lrz)lA  A1A 
^rz^lB  A1B  " ^OO^l  A1 


^zz’lA  A1A  ' (azz}  IB  A1B  + (Trz)  lR  AlR  " ^rz'lL  AlL  >(6) 


F3  ^rr^R  A2R  “ (°rr)2L  A2L  + (lrz}  2A  A2A 

|Trz’ 2B  A2B  ~ (O00) 2 A2 


F4  " zz)  2A  A2A  " (ozz}2B  A2B  + {xrz}2R  A2R  " (Trz*2L  A2L 


where  all  quantities  to  the  right  of  the  equal  signs  are  eval- 
uated at  time  tn , and  the  A's  denote  partial  areas. 

Upon  requiring  equality  between  the  normal  components 
of  velocity  for  both  materials,  one  obtains  the  expression: 

uf1  sin0  - Vj+1  cos0  = sinO  + V^+1  cos0  = 0 (7) 

To  complete  the  above  system  of  equations,  the  dependence 
of  the  frictional  shear  stress  x on  the  other  independent  var- 
iables must  be  prescribed.  In  a related  DNA-sponsored  program, 
a series  of  experiments  was  performed  by  S3  to  characterize  the 
frictional  shear  stress  developed  between  steel  and  rock  during 
high  speed  sliding . ^ Experimental  results  were  obtained  in 
this  study  for  sliding  velocities  up  to  3000  cm/sec  and  for 
normal  stresses  up  to  about  1 kbar . Within  this  range  of  slid- 
ing velocities  and  normal  stresses,  the  frictional  shear  stress, 
t,  in  all  cases  studied  was  found  to  obey  the  law: 

t = V Jn  exp  I - ~J2-  J ( 8 ) 

where  C denotes  the  sliding  velocity,  a is  the  stress  normal 
to  the  direction  of  sliding,  and  u and  t,  are  mater ial -dependent 
parameters.  For  the  class  of  defense-related  impact  problems 
of  interest  to  the  present  study,  the  sliding  velocities  and 
normal  stresses  lie  well  beyond  the  range  covered  by  the  ex- 
periments noted  above.  In  such  cases,  it  can  be  argued 
that  the  frictional  resistance  will  most  likely  be  less  than 
that  indicated  by  the  low  velocity-low  stress  laboratory  tests, 
and  also  exhibit  only  a weak  dependence  on  sliding  velocity. 

On  this  basis,  the  following  expression  was  proposed  in  Refer- 
ence 11  and  adopted  in  the  present  study  to  describe  the  fric- 
tional shear  stress  t developed  between  steel  and  Dakota 
sandstone  for  high  velocity-high  stress  conditions: 


t = a on  exp  (-cn/b)  for  on>0  C9) 


Here,  a = 0.19  and  b = 1,6  kbars.  Equation  (9)  prevents  an 
unlimited  buildup  of  frictional  stress,  reaching  a peak  value 
of  0.112  kbars  at  on  = 1.6  kbars.  In  the  CAPRI  numerical 
scheme,  Equation  (9)  is  described  by  the  following  expression: 


t = sign  { (V") 2 - (Vfc) l } F(on) 


which  provides  the  proper  siqn  of  t in  terms  of  the  tangential 

velocities.  The  function  F(o  ) is  defined  by: 

n 


F(a  ) = aa  exp  (-a  / b) 
n n n 


where  the  values  of  the  constants  a and  b are  listed  above. 

Equations  (2)  to  (5) , (7)  and  (10)  form  a system  of  six 

linear  equations  which  is  solved  in  the  CAPRI  code  by  matrix 
inversion  to  yield  the  values  of  the  six  unknowns  + ^ , V™4 1 , 

U ^ ^ °n  an<^  T at  t^ie  t,nc5  t^ie  bagrangian  phase  of  the 

cal cul at  ion . 


2.1.2  Calculation  of  the  Deviatoric  Strain  Rates  for 

M it>  r : a 1 s i r.  S 1 ; p ' '■  '1 1 s 

As  noted  above,  the  preceding  equations  apply  to  the 
Lagrangian  phase  of  the  Eulerian  numerical  scheme , and  they 
ensure  that  the  basic  conditions  which  apply  at  material  inter- 
faces [see  Eq.  (1)J  are  satisfied  during  this  phase  of  the  cal- 
culation. In  order  to  equilibrate  the  normal  stress  components 
after  the  mass  transport  phase,  one  must  calculate  the  strain 
rates  for  each  material  in  the  slip  cell.  This  is  not  a trivial 
task,  since  the  material  interface  is  a surface  across  which 
the  strain  tensor,  and  consequently  the  deviatoric  strain  rates, 
are  discontinuous. 

A procedure  was  developed  for  the  CAPRI  code  which  pro- 
vides an  approximate  method  for  calculating  the  deviatoric 
strain  rates  in  each  slip  cell  material;  this  approach  utilizes 


the  velocity  components  of  a material  element  in  the  slip  cell, 
as  well  as  those  of  the  same  material  in  neighboring  cells,  to 
calculate  the  deviatoric  strain  rates  for  that  material.  The 
procedure  is  clearly  most  accurate  when  the  material  interface 
passes  through,  or  near,  the  cell  center,  and  becomes  more 
approximate  as  the  position  of  the  interface  deviates  further 
from  this  condition. 

2.1.3  Egu i 1 ibrat ion  of  Normal  Stress  Components  After 

Mass  Transport  pfrase 

After  completion  of  the  mass  transport  phase,  the  code 
utilizes  an  iteration  procedure  to  equilibrate  the  normal  com- 
ponents of  stress  in  the  two  slip-cell  materials.  In  the  iter- 
ation scheme,  the  deviatoric  strain  rates,  determined  by  the 
method  described  in  Section  2.1.2,  remain  fixed  and  the  den- 
sities of  the  two  materials  in  the  mixed  cell  are  adjusted 
until  the  components  of  stress  normal  to  the  interface  are 
equal.  The  present  iteration  scheme  is  particularly  tailored 
to  treat  the  case  of  interest  herein  in  which  one  of  the  slid- 
ing materials  is  a metal  and  the  other  a rock  material.  A 
rock-type  constitutive  model  introduces  considerable  complex- 
ities into  the  iteration  procedure  because,  for  such  materials, 
the  hydrostatic  and  deviatoric  components  of  stress  are  usually 
coupled  through  the  dependence  of  the  yield  strength  on  pres- 
sure . 

After  the  procedure  for  treating  frictional  sliding  be- 
tween materials,  described  in  the  preceding  sections,  had  been 
incorporated  into  the  CAPRI  code,  a simple  two-material  prob- 
lem, for  which  a theoretical  solution  could  be  obtained  to 
check  the  numerical  results,  was  analyzed  with  the  new  code. 

A description  of  this  analysis,  and  the  corresponding  results 
obtained  are  provided  in  Appendix  B.  As  shown  there,  the 
agreement  between  the  theoretical  and  numerical  solutions  is 
excellent . 
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2.2 


MASS  TRANSPORT 


2.2.1  Shortcomings  of  Existing  Methods 

( 9 ) 

In  state-of-the-art  Eulerian  codes  such  as  HELP, 
the  mass  flux,  of  material  i through  a boundary  of  a 

mixed  material  cell  is  usually  given  by  an  expression  of  the 
form: 

Am.  = p . v.  A . A t (12) 

1111 

where  denotes  the  donor  cell  density,  v^  is  an  interpolated 
velocity  representing  the  velocity  normal  to  that  boundary  at 
the  end  of  the  time  step  At,  and  A ^ represents  that  portion  of 
the  area  of  the  cell  boundary  intersected  by  material  i.  The 
area  A^  is  determined  by  tracer  particle  positions. 

There  are  several  shortcomings  in  the  mass  transport 
algorithm  described  above.  In  the  first  place,  substantial 
inaccuracies  can  arise  in  the  calculation  of  mass  fluxes  as 
a material  interface  crosses  over  an  entire  cell  boundary 
during  a single  time  step.  To  illustrate  this,  consider 
Figure  4 which  shews  the  positions  of  a material  interface 
before  and  after  it  crossed  a cell  boundary  at  times  t and 
t+At,  respectively.  According  to  the  interface  position  at 
time  t,  the  cell  boundary  is  intersected  only  by  material  2. 
Thus,  on  the  basis  of  the  algorithm  discussed  above,  the  area 
for  material  1,  A. , is  zero  and  no  material  1 is  transported 
across  the  boundary.  It  is  obvious,  however,  from  the  position 
of  the  interface  at  time  t+At  that  some  material  1 has  entered 
the  cell.  Moreover,  it  is  necessary  to  reduce  the  flux  of 
material  2 as  calculated  by  Eq.  (12)  since,  for  this  case, 
there  is  generally  less  mass  of  material  2 in  the  cell  than 
is  calculated. 

Another  difficulty  may  occur  if  an  interface  leaves  a 
cell  and  the  flux  of  material  calculated  by  Eq.  (12)  is  in- 
sufficient to  evacuate  all  of  that  material  from  that  cell, 
i.e.,  some  of  the  material  is  left  behind.  In  this  case,  specia 
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procedures  are  usually  employed  to  increase  the  flux  so  that 
all  of  the  material  is  evacuated. 

Since  the  difficulties  noted  above  vanish  as  the  time 
step  becomes  very  small,  subcycling  procedures  have  been 
introduced  to  reduce  errors  in  transport  across  mixed-material 
cell  boundaries.  This  is  accomplished  by  using  N time  steps 
for  mass  transport  calculation  within  each  normal  time  step 
At.  The  number  N is  selected  by  the  user,  and,  in  practical 
applications,  is  usually  taken  to  fall  in  the  range  2 <_  N <.4. 
One  disadvantage  of  subcycling,  of  course,  is  that  it  can 
increase  significantly  the  cost  of  a calculation.  Furthermore, 
although  the  errors  in  transport  are  substantially  reduced, 
they  are  not  eliminated;  they  are  usually  small  and  are 
quickly  damped  by  the  spurious  numerical  diffusion  character- 
istic of  Eulerian  codes. 

2.2.2  New  Algorithm  for  Mass  Transport 

In  the  course  of  the  present  program,  it  was  found 
that  the  mass  transport  algorithm  described  above  was  unsatis- 
factory for  the  class  of  impact  problems  of  interest.  In  par- 
ticular, the  relatively  low  stresses  generated  by  projectile- 
rock  impacts  were  significantly  perturbed  by  the  relatively 
large  numerical  errors  generated  as  material  interfaces  crossed 
cell  boundaries.  Futhermore,  the  history-dependent  parameters 
in  the  sandstone  model  made  it  "remember"  the  non-physical 
states  it  passed  through,  thus  further  degrading  the  reliabil- 
ity of  the  results  obtained.  Since  the  number  of  "subcycles" 
of  the  multimaterial  cell  mass  transport  algorithm  necessary 
to  reduce  these  errors  to  an  acceptable  level  would  have  made 
the  calculation  prohibitively  expensive  to  do,  it  was  decided 
to  attempt  to  improve  the  algorithm. 

The  improvement  essentially  involves  a redefinition  of 
the  area  term,  An,  in  Equation  (12)  whenever  either  of  the  two 
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In  each 


problem  cases  discussed  previously  are  encountered, 
case,  Am^,  the  mass  flux,  is  considered  to  be  known,  i.e., 
it  is  the  total  amount  of  mass  of  material  n remaining  in  the 
cell.  Then  the  area  term  is  defined  to  be  the  value  necessary 
to  satisfy  Equation  (12): 

A = Am  /p  v At . (13) 

n n n n 

The  area  terms  for  any  other  materials  which  intersect  that 
cell  boundary  are  now  proportionally  adjusted  to  satisfy  the 
constraint : 


EA  = A (14) 

n 
n 

where  A is  the  total  area  of  the  cell  face. 

The  algorithm  described  above  was  incorporated  into 
the  CAPRI  code  a'nd  used  along  with  "subcycling"  on  several 
test  problems  involving  vertical,  horizontal  and  diagonal 
movement  of  a solid  bar  in  uniform  motion  surrounded  by  a 
free  surface.  Virtually  no  numerical  noise  was  generated 
by  the  improved  mass  transport  algorithm. 

2.3  MATERIAL  FAILURE 

The  problems  considered  in  the  present  stridy  are  char- 
acterized by  extensive  failure  of  material.  The  penetration 
of  a projectile  into  rock  is  accomplished  through  extensive 
fracturing  and  comminution  of  the  rock.  Furthermore,  the 
severe  deformations  that  soft  nose  caps  will  likely  experi- 
ence during  penetration  can  be  expected  to  far  exceed  the 
failure  threshold  of  the  material. 

Regions  of  failed  material  move  through  an  Eulerian 
grid  as  a result  of  (a)  further  cracking  and  failure,  and 
(b)  advection  of  material.  An  approximate  method  was  devel- 
oped and  incorporated  into  the  CAPRI  code  to  (a)  reduce  the 
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effect  of  numerical  diffusion  on  the  identification  of  regions 
of  failure  and  (b)  allow  for  the  effect  of  failure  on  material 
properties.  This  was  accomplished  through  the  introduction  of 
a flow  field  fracture  variable  F.  This  variable  is  a material 
property  which  characterizes  the  failure  state  of  a material. 
If  the  material  in  a zone  is  unfractured,  F = 1,  while  F = 0 
for  a zone  occupied  by  fractured  material  only.  For  zones 
containing  both  fractured  and  unfractured  material,  the  value 
of  F lies  somewhere  between  zero  and  unity.  Since  F is  a 
material  property,  it  is  transported  through  the  Eulerian  grid 
in  the  same  manner  as  the  stress  deviators. 

As  presently  structured,  the  CAPRI  code  assumes  that 
failed  material  is  unable  to  support  deviatoric  stresses. 
Accordingly,  in  zones  containing  both  failed  and  unfailed 
material,  i.e.,  where  0 i F i 1,  the  shear  modulus  G and  yield 
strength  Y are  reduced  proportionately  in  the  following  manner 
to  reflect  the  degradation  in  material  strength: 

G = FG„ 


Y = FY 

0 

Here,  F denotes  the  fraction  of  failed  material  in  the  zone, 
and  Gq  and  Yq  represent,  respectively,  the  initial  unfailed 
values  of  the  shear  modulus  and  yield  strength. 


NUMERICAL  STUDY  OF  NOSE-CAPPED  PROJECTILE  IMPACTING  ROCK 


3. 


A two-dimensional  numerical  study  was  performed  with  the 
new  CAPRI  code  to  examine,  and  to  demonstrate,  its  ability  to 
analyze  the  response  of  highly  deformable  projectiles  striking 
rock.  The  details  of  the  study  are  given  in  this  section. 

3.1  DESCRIPTION  OF  IMPACT  PROBLEM 

Figure  5 depicts  the  initial  axisymmetric  configuration 
of  the  pro jectile-rock  impact  problem  considered  in  the  numeri- 
cal study.  The  projectile  consisted  of  a solid  cylinder  of 
D6A-C  steel  to  which  a soft  magnesium  alloy  nose  cap  of  the 
same  diameter  was  attached.  The  diameter  of  the  projectile 
was  2 cm  and  its  overall  length  was  4 cm.  The  nose  cap  was 
1 cm  thick.  It  was  assumed  in  the  numerical  study  that  the 
nose  cap  and  the  main  body  remained  permanently  bonded  together 
during  the  penetration  process;  frictional  sliding  took  place 
only  between  the  sandstone  and  projectile.  The  impact  velocity 

4 

was  6 x 10  cm/sec.  The  numbered  locations  in  Figure  5 denote 
the  initial  positions  of  Lagrangian  points  which  identify  the 
material  positions  at  which  time  histories  of  the  flow  field 
variables  were  monitored  throughout  the  calculation. 

The  Eulerian  grid  used  to  obtain  the  numerical  results 
reported  herein  is  shown  in  Figure  6.  In  all,  a total  of 
3,000  computational  zones  was  employed. 

3.2  CONSTITUTIVE  MODELS 

3.2.1  Projectile  and  Nose  Cap 

The  constitutive  behavior  of  the  projectile  and  nose  cap 
materials  was  described  by  simple  elastic-ideally  plastic  models 
Material  property  data  for  D6A-C  steel  was  obtained  from  Refer- 
ence 12  and  corresponding  data  for  the  magnesium  alloy  AZ80A-T5 
was  found  in  Reference  13.  Because  of  the  relatively  low 
stress  levels  expected  in  the  class  of  problems  of  interest, 
a simple  equation  of  state  of  the  form: 
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Figure  5. 


Initial  pro jectile-target  conf iquration 
showinq  the  gauqe  positions  where  stress 
histories  were  monitored.  The  dashed  line 
indicates  the  initial  position  of  the  inter 
face  between  the  projectile  body  and  the 
maqnesium  nose  cap. 
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was  employed.  Here,  K denotes  the  bulk  modulus,  p is  the  den- 
sity, and  pq  is  the  initial  value  of  the  density.  If  the 
accumulated  generalized  plastic  strain  in  a material  element 
reached  a prescribed  level,  the  material  was  assumed  to  frac- 
ture; when  this  occurred,  all  components  of  the  deviatoric 
stress  tensor  were  reduced  to  zero  (see  Section  2.3) . 

The  values  of  the  material  parameters  used  for  projectile 
materials  in  the  numerical  study  are  summarized  in  Table  1. 

3.2.2  Rock 


Dakota  sandstone  was  selected  for  the  numerical  study 
because  it  typifies  the  kinds  of  rock  materials  found  at  many 
hard  sites  of  defense  interest.  Furthermore,  the  complex  be- 
havior of  Dakota  sandstone  has  been  rather  extensively  investi- 
gated in  the  laboratory  by  the  U.S.  Army  Waterways  Experiment 

( 14 ) 

Station  and  Terra  Tek , and  recently  a constitutive  model 

for  this  material  was  developed  by  C RT . ^ ^ The  CRT  model 
describes  the  complex  plastic  behavior  of  Dakota  sandstone 
within  an  elastic-plastic  framework;  it  requires  knowledge 
of  a number  of  history-dependent  parameters,  including  the 
pressure,  the  extent  of  fracture  and  the  generalized  plastic 
strain.  Provision  is  made  in  the  model  for  describing  post- 
fracture response,  wherein  the  strength  of  the  fractured 
material  is  degraded  to  that  characteristic  of  the  fully 
crushed  condition.  In  most  regions  of  behavior,  the  model 
uses  an  associated  flow  rule;  however,  under  certain  conditions, 
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Table  1 


Material  Properties  of  Projectile  Materials 


Material 

Property 

Magnesium  Alloy 
AZ80A-T5 

j D6A-C  Steel 

Mass  density 
p (gm/cm3) 

1.8 

7.85 

Bulk  modulus 

2 

K (dynes/cm  ) 

3.89  x 1011 

1.6G3  x 10 

Shear  modulus 

2 

G (dynes/cm  ) 

1.71  x 101  1 

7 . 8 8 8 x 1 0 

Yield  stress 

2 * 

Yq  (dynes/cm  ) 

2. GO  x 10° 

13.0  x 101 

Generalized  plas- 
tic strain  frac- 
ture threshold. 

0.0142 

O 

o 

kD 

O 

* Y denotes  the  yield  stress  in  simple  tension. 


are  associated  with  given  (identifiable)  material  elements. 

In  a Lagrangian  formulation,  material  elements  maintain  their 
identity  throughout  a calculation,  and  the  task  of  keeping 
track  of  history-dependent  material  parameters  associated  with 
such  elements  poses  no  special  problems.  In  an  Eulerian  for- 
mulation, however,  the  task  is  more  complex  since  the  basic 
formulation  is  not  designed  to  maintain  the  identity  of  indiv- 
idual material  elements.  Because  of  this,  there  was  some 
uncertainty  at  the  outset  of  the  study  as  to  the  accuracy  that 
could  be  achieved  in  treating  such  a complex,  history-dependent 
material  model  with  an  Eulerian  framework,  and  a part  of  the 
effort  was  devoted  to  addressing  this  question. 

In  the  course  of  the  present  study,  we  found  that  the 
complex  CRT  sandstone  model  could  be  adequately  treated  in  the 
Eulerian  CAPRI  code  framework  if  the  six  history-dependent 
material  parameters  were  incrementally  updated  in  the  numerical 
scheme  by  the  same  approach  used  to  update  the  deviatoric 
stresses  (see  Ref.  9 for  details  of  procedure  for  updating 
deviatoric  stress  components) . To  demonstrate  the  accuracy  of 
this  procedure,  a one-dimensional  wave  propagation  problem  was 
formulated  as  a test  case,  and  analyzed  independently  by 

CRT,  using  their  WAVE-L  Lagrangian  code,  and  by  S3  , with 

the  Eulerian  CAPRI  code.  The  specific  details  of  the  impact 
configuration  considered  by  both  CRT  and  S3  are  given  in 
Appendix  A,  together  with  the  numerical  results  from  both 
calculations.  An  inspection  of  these  results  shows  that  the 
calculated  stress  wave  profiles  from  both  computations  are  in 
reasonable  agreement,  and,  thus,  demonstrates  the  accuracy 
with  which  the  Eulerian  CAPRI  code  can  treat  a complex, 
history-dependent  constitutive  model, 

3.3  EARLY  TIME  ONE-DIMENSIONAL  SHOCK  STATES 

Under  the  assumption  of  uniaxial  stress,  the  shocked 
states  produced  in  the  sandstone,  magnesium  and  steel  by  the 
impact  can  be  determined  analytically,  and  used  to  check, 


where  appropriate,  the  calculational  results.  A plot  of  the 
Hugoniots  for  sandstone,  magnesium  and  steel  appropriate  for 
the  material  configuration  and  impact  velocity  of  interest  is 
given  in  Figure  7.  As  this  figure  reveals,  the  sandstone 
Hugoniot  intersects  the  reflected  steel  Hugoniot  at  a stress 
of  approximately  26  kbars.  In  addition,  the  figure  shows  that 
the  shock  transmitted  from  the  sandstone-magnesium  interface 
produces  a shock  in  the  steel  of  approximately  29  kbars.  Thus, 
on  the  basis  of  one-dimensional  shock  propagation  theory,  which 
is  valid  only  at  early  times  near  the  axis  at  the  sandstone- 
magnesium  interface,  it  is  concluded  that  the  use  of  a soft 
nose  cap  causes  the  shock  transmitted  back  into  the  steel  main 
body  to  be  greater  than  if  there  were  no  nose  cap.  Conversely, 
if  the  nose  cap  had  been  fabricated  from  a material  which  is 
harder  than  steel,  such  as  a tungsten  alloy  for  example,  sim- 
ilar one-dimensional  considerations  indicate  that  the  stress 
propagated  into  the  steel  body  would  be  loss  than  if  (a)  the 
nose  cap  was  made  from  a material  softer  than  steel , and  (b) 
there  was  no  nose  cap.  This  result  follows  from  Figure  8, 
where  the  appropriate  Hugoniots  for  sandstone,  steel  and  a 
tungsten  alloy  are  shown. 

It  should  be  emphasized  that  the  above  conclusions  apply 
only  to  sandstone  targets,  and  to  other  rock  materials  that 
exhibit  Hugoniots  which  are  less  stiff  than  sandstone.  For 
rocks  having  sufficiently  stiffer  Hugoniots  than  sandstone, 
it  can  be  shown,  using  similar  arguments,  that  a soft  nose  cap 
serves  to  reduce  the  transmitted  stress  levels. 

Whether  or  not  the  above  conclusions  also  hold  when  two- 
dimensional  effects  are  taken  into  consideration  is  not  known; 
their  validity  could  be  investigated  by  performing  a similar 
numerical  calculation  with  the  CAPRI  code  for  a steel  projec- 
tile without  a nose  cap,  and  comparing  the  stress  levels 
generated  in  the  steel  with  those  determined  for  the  magnesium 
capped  projectile  in  the  present  study. 
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NUMERICAL  RESULTS 


4 . 

The  results  obtained  from  the  numerical  study  described 
in  Section  3 are  presented  here.  As  noted  earlier,  the  numer- 
ical calculation  was  taken  out  to  27  usee  after  impact.  At 
this  time,  significant  deceleration  of  the  projectile  had  not 
taken  place,  but  the  soft  magnesium  nose  cap  had  undergone 
such  severe  plastic  flow  that  its  thickness  had  become  smaller 
than  the  size  of  the  calculational  Eulerian  zone;  this  resulted 
in  three  materials  being  present  in  a single  zone,  and  in  the 
current  state  of  development,  CAPRI,  as  well  as  Eulerian  codes, 
is  unable  to  adequately  cope  with  such  a situation.  The  numer- 
ical calculation,  however,  could  have  been  continued  by  re- 
zoning the  Eulerian  grid  with  smaller  zones  and  continuina  the 
calculation,  but  this  was  not  attempted  due  to  limitations  of 
remaining  funds. 

4.1  CONFIGURATION  OF  SYSTEM  AT  VARIOUS  TIMES  AFTER  IMPACT 
Figures  9 and  10  depict  the  calculated  configurations  of 

the  projectile-rock  system  at  various  times  after  impact  up  to 
the  time  (27  usee)  at  which  the  calculation  was  discontinued. 

Also  shown  in  these  figures  is  the  progressive  development  of 
failure  regions  in  both  the  projectile  and  rock  materials. 
Although  significant  deceleration  of  the  projectile  had  not 
occurred  by  27  usee,  severe  distortion  and  flow  of  the  soft 
magnesium  nose  cap  has  taken  place.  By  3 y sec,  the  nose  cap 
had  completely  failed  and  its  thickness  was  reduced  to  a thin 
layer  of  flowing  metal  by  27  usee.  Note  also  that  some  failure 
has  also  taken  place  in  the  steel  at  this  time. 

4.2  ACCELERATION,  VELOCITY  AND  DISPLACEMENT  AT  SELECTED 
LAGRANGIAN  POINTS 

The  calculated  time  histories  of  the  axial  and  radial  com- 
ponents of  acceleration,  velocity  and  displacement  at  the  most 
interesting  Lagrangian  gage  points  noted  in  Figure  5 are  depicted 


in  Figures  11  to  22.  The  response  at  Gage  5 is  of  particular 
interest,  since  the  early-time  deformation  at  this  location  is 
uniaxial  strain  until  relief  waves  arrive  from  the  free  sur- 
face. Comparison  of  the  calculated  peak  velocity  at  early  times 
shown  in  Figure  16,  with  the  analytically  determined  velocity 
given  in  Section  3.3  shows  good  agreement.  The  reason  for  the 
gradual  increase  in  the  axial  velocity  at  Gage  5 at  later  times 
(say  past  10  usee)  is  not  known,  however;  this  could  be  the 
result  of  complex  wave  interactions  within  the  magnesium,  or  be 
due  to  features  of  the  numerical  scheme.  The  accuracy  with 
which  the  code  can  treat  the  response  of  the  magnesium  nose  cap 
undoubtedly  decreases  as  the  magnesium  is  further  deformed 
into  a thin  layer.  Inspection  of  the  velocity  and  displacement 
histories  in  the  magnesium  at  Gage  5,  shown  in  Figures  16  and  17 
indicates  no  significant  deceleration  of  the  projectile  has 
occurred  by  27  usee.  However,  the  velocity  history  for  Gage  3, 
shown  in  Figure  12 , reveals  that  there  has  been  some  small  de- 
celeration of  the  steel  main  body. 

4.3  STRESS  HISTORIES  AT  SELECTED  LAGRANGI AN  POINTS 

The  calculated  time  histories  of  the  axial  stress  a , the 

radial  stress  a , the  hoop  stress  a,  and  the  shear  stress  t 

r v rz 

at  the  most  interesting  Lagrangian  points  are  presented  in 

Figures  23  to  28.  For  the  Lagrangian  points  which  lie  on  the 

axis,  the  conditions  a =a a and  x =0  hold.  Here,  and  in  that 

r 9 rz 

which  follows,  compressive  stresses  are  taken  as  positive. 

Consider  now  the  calculated  stress  histories  in  the  steel 
main  body  at  Gage  3,  which  are  shown  in  Figure  23.  The  early- 
time  axial  stress  response  shows  a peak  of  about  8 kbars  which 
is  considerably  smaller  than  the  29  kbar  stress  level  predicted 
in  Section  3.3  on  the  assumption  of  one-dimensional  deformation. 
It  is  thus  concluded  that  the  early  time  response  is  not  one- 
dimensional, but  is  influenced  strongly  by  the  arrival  of  relief 


waves.  At  later  times,  the  axial  stress  increases  up  to  about 
14  kbars  possibly  as  a result  of  complex  interaction  between 
waves  contained  within  the  steel  main  body. 

The  stress  histories  in  the  magnesium  at  Gage  5 are  shown 
in  Figure  25.  At  early  times,  the  response  at  this  location  is 
one-dimensional.  The  calculated  peak  axial  stress  when  one- 
dimensional conditions  prevail  is  15.8  kbars  which  compares  well 
with  the  value  of  17  kbars  determined  from  Figure  7 in  Section 
3.3.  The  magnesium  in  the  neighborhood  of  the  gage  fails  at 
about  1.5  ysec,  and  thereafter  behaves  as  a strengthless  solid. 
After  the  initial  pulse,  the  axial  stress  increases  up  to  a 
nominal  level  of  about  18  kbars  at  late  times.  The  cause  of 
this  rise  in  axial  stress  is  not  currently  known.  As  noted 
earlier  in  Section  4.2,  it  may  result  from  wave  interactions 
within  the  magnesium,  or  be  due  to  possible  deficiencies  in 
the  numerical  scheme  for  dealing  with  the  nose  cap  material 
as  it  deforms  into  a layer  which  has  a thickness  on  the  order 
of  the  zone  size. 

In  Figure  27,  the  stress  histories  in  the  sandstone  at 
Gage  7 are  shown.  Again,  the  calculated  peak  axial  stress  of 
about  8 kbars  at  early  times  is  substantially  smaller  than  that 
predicted  by  one-dimensional  considerations  in  Section  3.3. 

Thus,  the  early  time  response  at  this  location  also  appears 
to  be  significantly  affected  by  relief  waves. 

Finally,  the  stress  histories  at  the  off-axis  Lagrangian 
locations  are  given  in  Figures  24,  26  and  28. 

4.4  STRESS  CONTOURS  IN  THE  PROJECTILE  AT  SEVERAL  TIMES 

AFTER  IMPACT 

In  order  to  depict  the  evolution  of  stress  distribution  in 
the  projectile  contours  showing  the  stress  fields  in  the  pro- 
jectile at  5 , 15  and  25  usee  after  impact  are  presented  in 
Figures  29  to  31.  In  each  figure,  contours  of  the  radial  stress 
o^,  the  hoop  stress  a^,  the  axial  stress  o^,  and  the  shear  stress 
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TriT  are  shown  separately.  The  contour  intervals  are  two  kbars, 
except  in  the  case  of  the  shear  stress  where  they  are  one  kbar. 


Pro jectile-rock  configuration  at  various 
impact,  showing  region  of  failed  material 


36 


s <>•<*•  j 


.0  2.S  S.O  7.S  10.0  1J.S  'SO  17. S 70.3 

-tre  inionrsECTfCSi 


Figure  12.  Axial  velocity  history  at  Gage  3 


5.  CONCLUSIONS 


The  study  reported  herein  should  be  regarded  as  an 
initial  attempt  to  develop  a numerical  method  suitable  for 
analyses  of  problems  involving  the  impact  of  nose-capped 
projectiles  with  rock.  Such  problems  are  characterized  by 
large  plastic  deformation,  extensive  material  failure,  fric- 
tional sliding  between  cracked  (comminuted)  rock  and  the 
projectiles  metal  surfaces,  and  systems  involving  multiple 
materials . 

In  view  of  the  complexities  noted  above,  neither  the 
standard  Eulerian  nor  the  standard  Lagrangian  approaches  are 
well-suited  for  numerical  analysis  of  such  problems.  A deci- 
sion was  made  at  the  beginning  of  the  effort  to  adopt  the 
Eulerian  approach  as  the  basis  on  which  to  develop  a numerical 
method  for  this  class  of  problem;  this  decision  was  difficult 
and  made  only  after  weighing  the  many  advantages  and  disadvan- 
tages of  both  approaches.  The  Eulerian  approach  appeared  to 
offer  more  promise  primarily  because  of  its  superior  ability 
to  deal  with  large  deformations,  which  we  felt  to  be  important 
in  the  analysis  of  highly  deformable  projectiles.  Lagrangian 
methods  are  nevertheless  frequently  applied  by  the  defense 
community  to  problems  involving  large  deformations, and  numer- 
ical solutions  are  made  possible  through  extensive  rezoning  proce- 
dures. Little  can  be  said,  however,  about  the  credibility 
(accuracy)  of  results  so  obtained. 

The  main  objectives  of  this  study  were  accomplished, 
namely  (a)  a new  Eulerian  code  CAPRI  was  developed  which  is 
specifically  tailored  to  treat  problems  involving  the  normal 
impact  of  deformable  projectiles  with  rock,  and  (b)  the  appli- 
cation of  CAPRI  to  a projectile-rock  impact  problem  was  demon- 


The  CAPRI  code  is  a two-dimensional  Eulerian  code  devel- 
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oped  from  the  existing  HELP  code  framework.  It  features  the 
following  new  capabilities  which  are  evidently  important  for 
the  class  of  impact  problems  of  interest  to  the  present  inves- 
tigation . 

• Accounts  for  material  strength  in  mixed- 
material  computational  cells. 

• Treats  frictional  sliding  between  material 
surfaces  in  contact. 

• Attempts  to  satisfy  the  correct  boundary 
conditions  on  stress  and  velocity  at  inter- 
faces between  materials  in  contact. 

• Contains  an  improved  algorithm  for  trans- 
porting mass  and  other  related  material 
properties  through  the  Eulerian  grid. 

• Accommodates  complex  geologic  material  models, 
including  those  having  a number  of  history- 
dependent  parameters  (see  Appendix  A for 
further  discussion  on  this  point) . 

The  major  portion  of  study  was  devoted  to  the  formulation, 
development  and  debugging  of  the  intricate  numerical  algorithms 
required  to  incorporate  the  features  listed  above. 

To  examine  and  demonstrate  the  ability  of  the  new  code 
CAPRI  to  perform  numerical  analysis  of  deformable  projectile- 
rock  impacts,  it  was  applied  to  a problem  involving  a solid 
cylindrical  projectile,  having  a soft  magnesium  nose  cape,  im- 
pacting sandstone  at  normal  incidence.  A numerical  solution, 
which  appears  to  be  acceptable,  was  obtained  out  to  a time  after 
impact  (27  ysec)  at  which  a substantial  part  of  the  penetration 
process  had  occurred.  In  the  process  of  performing  this  calcu- 
lation, however,  a number  of  difficulties  with  the  numerical 
scheme  arose  which  made  it  impossible  to  continue  the  calcu- 
lation until  the  source  of  the  difficulty  was  uncovered.  An 
appropriate  modification  was  made  in  the  flow  field  at  that  time, 
and  the  calculation  was  resumed.  As  noted  earlier,  however,  the 
modifications  were  of  limited  extent  and  made  in  material  regions 
where  their  effects  were  probably  of  negligible  significance. 


The  calculation  was  not  continued  beyond  27  ysec  because, 
by  this  time,  the  soft  magnesium  nose  cap  material  had  under- 
gone such  severe  plastic  flow  that  its  thickness,  in  some 
region  became  less  than  the  Eulerian  zone  size;  this  resulted 
in  three  materials  being  present  in  one  zone  and,  in  the  present 
state  of  development,  CAPRI  is  unable  to  adequately  cope  with 
such  a situation.  The  calculation  could  have  been  continued 
further  by  rezoning  the  calculational  field  with  smaller  Euler- 
ian zones;  this,  however,  would  provide  only  a temporary  so- 
lution if  further  thinning  of  the  magnesium  layer  reduced  its 
thickness  to  that  of  the  new  zone  size.  Two  other  alternatives 
appear  to  be  more  desirable,  namely,  either  generalize  the 
numerical  scheme  to  handle  (admittedly  approximate)  three 
materials  in  a cell,  or  in  very  thin  magnesium  regions  remove 
the  magnesium  and  reposition  the  sandstone  interface  to  be  in 
immediate  contact  with  the  steel. 


Some  general  comments  on  the  study  are  as  follows; 

(1)  No  firm  conclusion  can  be  drawn  from  the  study 
reported  here  regarding  the  effectiveness  of  nose  caps  in  re- 
ducing the  level  of  stress  propagated  back  into  the  main  body 
of  a projectile  impacting  rock.  To  address  this  question, 
numerical  calculations  should  be  performed  with  the  CAPRI  code 
for  projectiles  having  (a)  no  nose  cap  and  (b)  nose  caps  of 
material  harder  than  the  steel  main  body,  say  a tungsten  alloy. 
On  the  basis  of  one-dimensional  response  considerations  pre- 
sented in  Section  3.3,  however,  it  was  found  that  the  use  of 
a soft  nose  cap  served  to  increase  the  level  of  stress  propa- 
gated into  the  main  steel  body.  Conversely,  the  use  of  a 
nose  cap  of  material  harder  than  steel  serves  to  reduce  the 
transmitted  stress  level.  The  validity  of  these  conclusions, 
based  on  one-dimensional  considerations,  needs  to  be  investi- 
gated under  more  general  two-dimensional  conditions. 
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(2)  The  numerical  approach  developed  and  incorporated 
in  CAPRI  for  treating  frictional  sliding  between  materials 
appears  to  be  reasonably  accurate  so  long  as  the  volumes  of 
the  two  sliding  materials  in  a mixed  cell  are  approximately 

of  the  same  order  of  magnitude.  However,  as  the  sliding  inter- 
face moves  through  the  fixed  Eulerian  grid  during  the  calcu- 
lation, situations  will  frequently  arise  wherein  the  volume 
of  one  material  is  very  much  less  than  that  of  the  other.  As 
the  relative  volume  of  one  of  these  materials  decreases,  the 
accuracy  of  the  present  numerical  algorithm  for  treating  the 
motion  of  the  smaller  volume  material  also  decreases.  This 
may  not,  however,  cause  a significant  problem  since  a numeri- 
cal error  made  in  a small  volume  of  material  would  be  expected 
to  have  much  less  effect  on  the  overall  response  of  the 
projectile-rock  system  than  a comparable  error  made  in  a 
larger  volume  of  material . 

(3)  It  is  not  clear  from  the  results  of  the  present 
study  that  the  additional  complexity  introduced  into  the  CAPRI 
code  by  the  incorpora tion  of  the  sliding  friction  capability 
is  justified,  since  the  levels  of  shear  stress  generated  by 
high  speed  sliding  between  steel  and  rock  are  evidently 
relatively  small.  The  frictional  sliding  law  adoDted  in 

the  calculation  presented  herein  was  suggested  to  us  by  Dr. 

Paul  Hadala , U.S.  Army  Waterways  Experiment  Station,  on  the 
basis  of  results  from  a related  experimental  study. ^ ^ 
According  to  this  law,  the  peak  shear  stress  that  can  be 
generated  from  frictional  sliding  is  about  0.1  kbars  and,  for 
the  range  of  normal  stresses  of  major  interest,  the  shear 
stress  is  expected  to  be  considerably  less  than  0.1  kbars. 

It  is  not  clear  what  influence  such  relatively  small  shear 
stresses  have  on  (a)  the  nose  cap  deformation  and  (b)  the 
overall  projectile  deceleration.  If  their  effect  is  small 
and  can  be  reasonably  neglected,  it  would  provide  considerable 
simplification  in  the  numerical  scheme;  this  was  not  investi- 
gated, however,  in  the  present  study. 


(4)  The  ability  of  the  CAPRI  code  to  treat  problems 
involving  the  impact  of  highly  deformable  projectiles  with 
rock  was  demonstrated  for  the  case  of  a projectile  that  had 
a solid  body.  Although  the  scope  of  this  study  did  not 
include  analysis  of  several  nose  cap  materials,  a great  deal 
can  be  learned  about  the  effectiveness  of  various  nose  cap 
materials  from  analysis  of  nose  capped  solid  projectiles. 

With  the  existing  version  of  the  CAPRI  code,  however,  calcu- 
lations can  be  performed  for  projectiles  having  shapes  of 
more  practical  interest.  Rather  than  the  blunt  nose  type 
projectile  considered  in  the  present  study,  numerical  studies 
should  be  performed  for  ogive-shaped  projectiles  having  a 
variety  of  nose  cap  materials  in  order  to  fully  understand 
the  role  that  nose  caps  play  in  effecting  stress  levels  trans- 
mitted into  the  main  body  of  a projectile.  Before  such  studies 
are  undertaken,  however,  we  recommend  that  some  effort  be 
devoted  to  improving  several  portions  of  the  code. 
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APPENDIX  A 


COMPARISON  BETWEEN  LAGRANGIAN  AND  EULERIAN  CODE 

CALCULATIONS  FOR  ONE-DIMENSIONAL  WAVE  PROPAGATION 

IN  SANDSTONE 

In  order  to  examine  how  well  the  complex  CRT  sandstone 
model,  with  its  six  history-dependent  parameters,  can  be  treated 
by  the  Eulerian  CAPRI  code,  a one-dimensional  test  problem 
involving  wave  propagation  in  sandstone  was  formulated and 
numerically  analyzed  with  both  a Lagrangian  code  and  an  Eulerian 
code.  The  problem  consisted  of  a sandstone  flyer  plate  impact- 
ing a target  of  like  material.  The  initial  configuration  of  the 

system  is  shown  in  Figure  32.  The  impact  velocity  was  chosen 
4 

to  be  3 x 10  cm/sec  so  that  the  impact  stresses  initially 
developed  in  the  flyer  and  target  would  be  about  8 kbars ; at 
this  stress  level  the  deviatoric  component  of  stress  is  a sub- 
stantial part  of  the  total  stress,  and  it  is  not  masked  by  a 
much  larger  hydrostatic  component. 

A uniform  zone  size  of  0.05  cm  was  used  in  both  the 
target  and  flyer  plate.  The  calculation  was  taken  out  to  2.5 
x 10  sec  to  permit  the  formation  and  subsequent  propagation 
of  the  tensile  wave  the  develops  near  the  rear  surface  of  the 
target;  this  permitted  the  code  to  exercise  the  tensile  cutoff 
feature  of  the  CRT  sandstone  model.  However,  no  separation  of 
the  material  was  permitted  in  the  calculation  when  the  tensile 
cutoff  was  invoked. 

At  S^'s  request,  CRT  performed  a numerical  analysis  of 
the  above  problem  with  their  Lagrangian  WAVE-L  code,  to  provide 
us  with  Lagrangian  results  against  which  we  could  compare  the 
corresponding  Eulerian  CAPRI  code  results.^6'  The  results  of 
the  CRT  numerical  study  are  shown  in  Figures  33  to  37  where  the 
calculated  stress  and  generalized  plastic  strain  fields  at  a 
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Figure  32.  Initial  configuration  of  one- 

dimensional  impact  test  problem. 
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number  of  selected  times  after  impact  are  given.  The  results 
obtained  by  using  the  new  Eulerian  CAPRI  code  are  depicted 
in  Figures  38  to  42,  where  calculated  stress  and  generalized 
plastic  strain  fields  at  the  same  times  as  the  CRT  results  are 
shown.  A comparison  of  the  corresonding  numerical  solutions 
reveals  that  there  is  good  agreement  between  the  stress  and 
generalized  plastic  strain  fields  predicted  by  the  two  differ- 
ent calculational  approaches,  and  thus  demonstrates  that  the 
spurious  numerical  diffusion  of  an  Eulerian  approach  is  not  a 
serious  factor,  at  least  at  these  stress  and  velocity  levels. 


predicted  by  CAPRI 


predicted  by  CAPRI  code. 
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NUMERICAL  ANALYSIS  OF  A SIMPLE  PROBLEM  INVOLVING 
FRICTIONAL  SLIDING  BETWEEN  TWO  SURFACES 

The  configuration  of  the  problem  considered  in  this  sec- 
tion is  depicted  in  Figure  43.  Here,  a cylindrical  rod,  in- 
itially at  rest,  is  put  in  motion  through  sliding  frictional 
contact  with  an  annular  region  of  like  material  moving  with  an 

initial  axial  velocity  of  V . As  a result  of  the  friction  con- 

o 

tact,  the  axial  velocity  of  the  rod,  V^ , increases  while  that 
of  the  outer  material,  V^ , decreases.  Eventually,  the  two  mater- 
ials reach  the  same  axial  velocity. 

The  materials  are  assumed  to  behave  elastically  with  the 
bulk  and  shear  moduli  taken  to  have  the  respective  values  of 
1,280  and  1,600  kbars . The  initial  density  of  the  material  was 
taken  to  be  7.8  gm/cm^ , while  the  frictional  shearing  stress 
to  between  the  sliding  materials  was  assumed  to  have  a constant 
value  of  1.22  kbars.  The  inner  material  was  initially  at  rest, 
while  the  outer  material  was  given  an  initial  velocity  of  1.0 
x 104  cm/sec. 

B. 1 THEORETICAL  SOLUTION 

A theoretical  analysis  of  the  motion  of  the  two  materials 


acted  upon  by  a constant  frictional  shear  stress,  x , at  their 
common  interface  leads  to  the  following  expressions  for  the  var- 
iation of  the  axial  velocities  V^  and  V?  with  time: 


Transmittive 
Boundary 
(both  ends) 


Interface  Between 
Materials  with 
Frictional  Sliding 


Computational  Grid 
Used  in  CAPRI  Code 
Calculation.  The 
Dimensions  of  all 
Zones  Were  Ar  = Az 


z = 1 cm 


Frictionless , 
Rigid  Boundary 


Figure  43.  Conf ig"ration  of  problem  considered  to  check  out 
new  sliding  interface  capability  in  CAPRI  code. 
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where  o is  the  density  of  the  two  materials , and  r^  and  r0 

are  defined  in  Figure  43.  Furthermore,  it  can  be  shown  that  the 

two  materials  reach  a common  velocity  of 


at  time,  t^ , given  by: 


B . 2 NUMERICAL  RESULTS 

In  the  numerical  solution  to  the  problem  depicted  in 
Figure  43  a uniform  Eulerian  grid  having  zones  with  dimensions 
ir  = iz  = 1 cm  was  used.  As  indicated  in  Figure  43,  three 
computational  zones  were  used  in  the  radial  direction,  with  the 
(sliding)  material  interface  passing  midway  through  the  interior 
zone.  In  the  axial  direction,  eight  zones  were  employed,  and 
transmittive  boundaries  were  used  to  introduce  and  remove  mat- 
erial from  both  ends. 

The  CAPRI  calculations  were  carried  out  to  a time  slightl 
less  than  that  at  which  the  two  materials  reached  a common  vel- 
ocity. The  results  of  the  calculation  are  depicted  in  Figure 
44,  where  the  variation  of  the  velocities  V,  and  V-,  with  time 
are  shown.  Also  shown  in  this  figure  are  the  variations  of 
V,  and  V ^ with  time  as  given  by  the  theoretical  expressions  in 
Equations  B.l  and  B.2.  As  an  inspection  of  Figure  44  reveals, 
the  agreement  between  the  CAPRI  code  numerical  results  and  the 
theoretical  solution  is  extremely  good  for  this  relatively  simpl 
problem  in  which  the  sliding  interface  passes  midway  through 
the  calculational  cells  and  remains  parallel  to  the  vertical 
grid  lines.  Similar  accuracy  cannot,  however,  be  expected  for 
more  complicated  problems  in  which  the  sliding  interface  is  not 
centrally  positioned  in  the  calculational  grid. 
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